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Abstract 


In the summers of 1963 and 1964 a double-ended 
refraction profile was shot between Suffield, Alberta and 
swift Current, Saskatchewan. A six layer earth model was 
obtained from these data, the crustal thickness under the 
east shot point (Swift Current) being 42 kms and the 
thickness under the west shot point (Suffield) being 51 kms. 
theseidata aressimilarsin general form to, those reported by 


MeyernandeMcCamy csineMontanas,+justseuthrofmour,lines 


Previous control west from Suffield has been ex- 
vendea inco the * Rocky Mountains. Since these data are not 
reversed, definitive interpretation is not possible. Some 


aluernative interpretations are presented. 


ppecural "analyses were also~-carrred out “on many of 
wie recorug., “rower spectra indicate thav~there*ts "little 
Chanse in Trequency with distance from the shot point. They 
do, however, indicate that there are two dominant frequencies 
eoLeseismiceinformation, onesateabout.5 cpsrsandysanother at 
apeurelOncpss. shobrenvinonment i\determines (whieh ofrthese 


two-is more important. 
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The Long Range Refraction Survey 


Ever since man has been inquisitive enough to 
SLuay Une nature of his surroundings, the structure of the 
earth has presented a major problem. His descriptive studies 
Oh whe surrace of the earth had revealed only superficial 
information about the structure of our planet. With the 
advent of the seismometer and the development of the science 
of seismology (a geophysical science which is concerned with 
the study of earthquakes and the measurement of the elastic 
properties of the earth), man has been able to gain a better 
Loe eoL bnvoe the) Structure of the interior of the earth. 
From these measurements of earthquakes the geophysicist 
Was Sable wwordistinguish three broad 'regisons) in thevearth, 
LHese Were: 
ihe) The truest ofthe earth.) This: Yegiom icomprises the oucer 


shell of the earth, and measures about 35 kms thick under 


Contiments Send Sbout. 20 to TPS kmetthicktander oceans, The 


base’ of the crustal layer is defined by ‘the Mohorovicic 
discontinuitye Thisitisrarselsmic:discontinuity,.iand is 
recognized throughout the world. 

De The Wantlee. This region of the earth extends down from 


thea pase on phe crust to.s depth of -900 kms. This layer 


4s believed to be composed of iron rich silicates, although 


near the Mohorovicic CLoconiUuaey 2G is vbelleved to be an 
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ultra-basic silicate much like dunite, peridotite or 

eclog tezsarcincesthetmantlerwill sustaintcboth «longitudinal 
andeeneariwaves sit iis: thought rho tbersolid, “ati least as 

far as seismic waves are concerned. 

Janebnceecore;, i cihbeweyohe sinvertorrof thesearth;- starting 
at uné centre and having a radius of 3470 kms. It is 
thought to be composed of iron and nickel in the free 
sGace, and furcnermore, is thovent to be In the fluid 
state, except for a small region in the centre known as 
tie Liner core, Tne reason for believing that It is 
ayy the fluid suave te the fact that 16 will not sustain 


shear waves in the seismic frequency region. 


Rarthquakes have provided a great deal of information 
AuGuuietne eartny but because: on the inability of the ‘seismolo- 
PUStwtO get sbnevexaeti time sand location“of occurrence, where 
has always been a certain degree of uncertainty in the data 
derived crezarding ttherrecrust of the*earth..\iMore recently, 
seismic equipment has become available with much higher 
resolution than the old seismometers used in earthquake 
measurements, and consequently controlled seismic experiments 


are now possible. 


Belsemic surveys of the long range type have only 
peen carried out since about. the early 1950's... The nature of 
these surveys has been to try and delineate the structure of 


Ene erust down LO the Mohorovicic discontinuity. It was expected 
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that the crust would be composed of different layers, having 
diilerent seismic properties, and as a consequence of this 
USsel Ul Serene retlecvion and retraction data could be ob- 
tained from these boundaries. . The implementation of long 


Vanse (Telractlon Surveys 18 reélatively straigntrorward. 


in, general, one. chooses. an, area,in, which. one wishes 
Dopvueverimilico. Loe, Crustal structure, .Afver selection of the 
eres, ono.  noles,are ori lied at both ends of, the, line, of 
ImLerest., and recording locatlons.are interspersed at 
uniform distances between the shot points. At prearranged 
times, explosive charges (usually TNT or dynamite) are 
exploded tuscue or Lie two. shot loles, the.charges providing 
a source of energy for the propagation of elastic waves 
through the crust. These waves are reflected and refracted 
along the various boundaries of the crust, and return to the 
Burtrace of the earth. At the recording locations, trucks 
equipped with instruments which are sensitive to ground 
motions (geophones) record these reflected and refracted 
arrivals. From the measured times of arrival of the waves, 
it is possible, using simple ray path methods, to determine 
Sne Structure and velocity distribution of the crust. see 


fie, A for a diagram illustrating this method. 


Such data are of fundamental importance to the 


seismologist. Long range surveys provide a good approximation 
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Or one slayered crustal structure, and where it is*®important, 
this layered structure may be interpreted in terms of the 
surface topography (e.g., under mountains one would expect 
UNerTe compe, some detrormation of the mantle at thesMoho 
boundary which would reflect the fact that mountains exist 
above it). Furthermore, these data allow us to make certain 
deductions about the elastic properties of the crustal 
materials, what the form of these materials might be, and 
whether the boundaries in the crust represent sharp discon- 
tinuities or simply transition zones, which are perhaps only 
a wavelength or so wide. Refraction surveys carried out on 
a local scale cannot hope to answer such broad questions, 
but with the ever increasing amount of data, a more universal 


picture should become available to the geophysicist. 
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Shot Hole 


Fig. A. 


Reflected Rays 
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Recordime Truck 


Geophones 


Refracted Rays 


Simple Diagram to Illustrate Ray Path Propagation 


in a- Layered Earth, and the Method of Recording. 
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The present, program of crustal studies at the 
URLVersity of Alberta had its inception in the, summer of 
1962. Since this time a number of experiments have been 
carried out and completed, among them being a double ended 
retraction profile between Suffield, Alberta and High River, 
Alberta, a double ended refraction profile between Suffield, 
Alberta and Swift Current, Saskatchewan, an extension of the 
SeuLtlelo data to"vne weet OuL to a distance “of ~500™kms., ‘arid 
the start of a program to detect near vertical’ reflections 
from tne déep crustal layers. Furthermore, “the University 
took part in the Lake Superior experiment which was conducted 
in the month of July, 1963. Weaver (1962) has analyzed 
tie waca ut rom Une .utitield - High River profile and reports 


Goie LOLLowLne results: 
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Depths under the East Shot Point (Suffield) 


to limestone (Miss. ) = tel, ne 
to Precambrian = 2.07 King 
to low velocity layer = 22,07. oie 
to Conrad = 35.25 nine 
to Moho = 48.16 kms 


the crue velocities tinweach ‘of these (layers swas as ‘follows: 


V, = 3.18 kms/sec Vy = 6.10 kms/sec (assumed) 
Vo = 5.88 kms/sec V, = 7-32 kms/sec 
V3 = 6.40 kms/sec Ve = 8.25 kms/sec 


Weaverswas forced to contrive a low velocity layer in order 
oO expla tuie 1act bDhat critical distance calculated from a 
Strole niuhorwardeinterpretacrion of the refraction.graphs 

Was much @reater than that.actually observed. He does, 
bowewerneoresent.a.2000 argument as to why there could 
Pause exis, Buch a:layer in southern Alberta; the reader 


is referred to Weaver (1962) for this argument. 


During the summers of 1963 and 1964 a double- 
ended refraction profile was completed between Suffield, 
Alberta and Swift Current, Saskatchewan. A total of 30 
records were obtained along this line, 17/7 from the east 
shot point (Swift Current) and 13 from the west shot point 


(Suffield). These records were plotted on reduced travel 
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tame @rapns,- and-inconjunction with three continuous 
velocity logs, aisix layersearthamodel waswpostulatedntor 
tne. Clusval='slructure in’ southern Saskatchewan. > Some prob= 
lems Gio vVarise inrthe anterpretationt of the! data... These 


Wate oS GPSCUSSEGLIN Cheilr appropriates section. 


Most. Of cine records alone bovh lines Were digitized, 
Pid. wie COWGCE SopeCL ra Were CcompuLled ror uhese records. 
fuWwee LOuUnG bile Bpec trum of the recordas was closely related 
to une geq0logic environment or the shot holes. Furthermore, 
Toei c Onin (piee. ay teas along the Surrleld cast line, 
there was little or no change in the frequency with 


distance. 


Peaveoues.on of the Switt Current = sufiield Line 
Well comperige the bulk or this thesis. There will be a 
Grier MeNtLon of the results of the westerly extension of 


ene pull iela -— Hien River line. 
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Chapter 1 


Ind) “WMstrumentation. 

pegorence tO figs. Ia aind 1b shows the instrumen- 
tatlom sysuem used at both the recording location and at the 
Siow oemvore= [he recording location systemiis relatively 
Simple. The signal from le different selismometers along the 
Spread are fed into twelve separate VLF refraction amplifiers. 
Cie 1BLGered; ampllincdssienal’ be.tnen fedeinto a recording 
Ceti tvorrapn, and a paperi record Groduceo eiethe signal. 
Aspone can see, there are three auxiliary signals fed into 
The -osciopranh. They are from a chronometer, a WWV. 
receiver, ang a standard broadcast receiver. These signals 
Were reproduced in order to facilitate the timing of the 


spreacg: Location records. 


Pie recording system at the shot, point 1s evenemore 
simple. Themshovueinstant from the shooting box, the signals 
from a chronometer, a WWV receiver, and a standard broadcast 
receiver are all fed into a recording oscillograph, and a 
peper record: producedsof ithese sigvais.|-This record™is Used 


oon Weekerhe shot instant onto! she, spreadcrecords. 


This system was used for most of the work which is 
dise@sseurin this thesis. |During the spring of 1964, however, 
the VLF amplifiers were modified in order to increase sensiti- 
vity and improve signal to noise ratio. At the same time, 
outputs were provided for tape recording. A few records were 


made using this modified system. 
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1.11 TIS-36 Seismometer. 
The Texas Instruments S-36 seismometer was used 
for most of the geophysical measurements carried out at the 


University of Alberta. Its characteristics are as follows: 


(a) The detector has a double coil system wound as a 
"hum bueking” device. 
(b) Tie ovira lPeereQieney 25 = coe’ = VO eps. 


(c) The coil resistance is 4000 ohms. 


ine Trequeney response characteristics are shown in 


aes 2 


ioe Vor herracvion Amplifier. 

Texas Instruments VLF refraction amplifiers have 
Ceriie“seo, Unroucniou, the operations. The total drain of 
bie eto VY power “supply is 55 milliamps. The system perfror- 


Meance, Dicures Tor the unmodified amplifiers are as Tollows: 


(a) The frequency response is 0.7 to 33 cps. Fig. 3 shows 
the system response curves for six different. filter 
settings. 

(bo) Fig. 4. gives the phase response curves in milliseconds. 

(ec) Distortion is below 1% at 7 cps, for inputs below 2 milli- 
volts, with the gain adjusted to give a 1" peak-to-peak 
deflection. Fig. 5 shows the variation of distortion with 


input under the conditions stated above. 
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POr ah inpua Of 15 microvolts, Zero db. attenuation and 

be ie SUD filter set at 46 eps, the deflection of the 

high gain trace is 1" peak-to-peak. See Fig. 6 for 
Overall magnitication of system, including geophones. 

The input ampedance is 10,000 ohms, and the line is 

balenced, No input transformer is used. 

Tne Internal 1 nstrumen, roles le equivalent to aninput 
fib niel Or s TM erovoice - 

TNers aré~dual™ourourse, -having a ‘derlectiron ratio of 
uta i, "NO oOULDUC tranerormer is-useda;:"one side of the 


GUGDUL LIne 1s” 2rounced., 


The external features of the system are as follows: 


(a) 


(b} 


Inere ore TWO 1nput. plugs, having precisely reversed 
polarities. 

OWeDChiNe Tact iIga vies anciude | line Tes. Switch vo eneck 
Tie COncibuipvy Of tne epreaa ling, separate avyrenustion 
switches for each amplifier (calibrated in 6 db. steps), 
Petit Ger secvite (SW..Cn Or Gach palr oF einplitiers 
(Calitrated With positions at 3, 12, 16, 24%, 32, and 

45 ops), and line balancing switches for each amplifier 
which siliipiy introduce resistive line balancing to create 
a common mode rejection of 60 cps signals induced by 


nearby power lines. 
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Oo. 
The modified amplifiers have the following added features: 
(a) The gain has been increased OM SS wacror of 15. 
(b) Input transformers are now used, maintaining an input 
Linpedance of 10,000 ohme. 
(ec) The line balancing circuit has been modified in its 


circuitry, but itcstilie performs» the samesoperation. 


there have been further radical changes. in the internal 
Mires try oO, one anmpliriers. These modliications were made 
in an effort to enhance the gain of the system and to 
enhance the signal to noise ratio. I1t appears that these 


mocirications have been eminently successful. 


wis Power Supply. 

The power supply is a fairly standard unit, 
Gediavyering TOOt volts DC. A le volt. de input is: inverved 
be m.c. DY 8 Transistorized inverter, anplitied, then 
Poecciried. ana filtered to give an output of 300 volts dc, 
This voltage is then passed through a tube type regulator, 
which finally delivers 100 + volts of well regulated dc. 
This voltage may be adjusted to suit the operating condition 


Sr cne amplifiers. 


This power supply was also modified to fit the new 
system of amplifiers. The output was broken up into three 
separate outputs. Regulation of two of the. outputs.is 
performed by the standard. voltage-regulator tube, circuit, 


while the critical reference voltage for the amplifiers is now 
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secu DynartO(i volt battery pack; rather than by V-R tube 
regulators. These voltages are not adjustable-by an exterior 


COLT ror. 


1.14 Recording Oscillograph. 
ine RS-6U recording oscillograph is a self 

contained unit composed of an optical system, a timing 

system, recording galvanometérs, and a paper transport 

Eyouvems. (1tUsS characteristics, are, 

a} Some 30 galvanometers having a natural period of 200 cps. 

(bo) Three paper speeds are available: SLOW inedium, wand fast, 
corresponding to speeds of about..2, 4,~and.10.inehes per 
second. 

(c) Timing lines are provided by a 100 cps transistor supply 
Senda Synchronous MoLor drivime .2-drumacontaining slits. 

(ad) There are 24 channels of seismic information possible, 
ber being the nich ‘gain: traces and -the.other 12 beine the 
low gain traces. Deflections are in the ratio of 4 tol. 
There are some six other auxiliary channels. 

(e) The intensity of the recording lamps is governed by the 
Paper -SceceadLswitch. Joel s orightes.£oreniehy speed, 
and cimmneso Tortchow speed. 

(f) he paner Uranspor..system-handles, a-roll of-sersmic 
photographic recording paper which is 8" x 200!. 

(2m Fixing, and developing of the records)is accomplished in 


the tanks fitted below the camera in a light tight box. 
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(h) Power requirements are 12 volts dc. Dissipation is at 
the rate of 13 amps during recording at the highest 
speed... The. timing.system draws,0.150.amps« «The RS -8U 


camera was not modified for the new system of amplifiers. 


ive-“PLeld° Procedure. 

The: trucks were.fitted with the following equip- 
ment<+-.,one.12) chagnel. TI VLE refractionfamplif terpyandg 
associated power supplies, a.recording oscillograph and icge 
associated developing tanks,-a,Westrex chronometer, -agWyv 
receiver, a standard broadcast receiver, 13 TI S-36 geophones, 
two geophone cable reels (each holding one mile of cable), 
assorted shovels..a tool. box,.spare- paper, and.additianad 


cans of developer and fixer. 


The field procedures adopted were fairly simple, 
calling for a maximization in time saving: and a®minimi sagi08 
in energy expenditure. Recording locations were usually 
chosen during planning sessions conducted in Edmonton. They 
were chosen with the following criteria in mind: first, 
Ghat each recording location ‘should 1ié46n, -or very oloese 
to the straight line which joined the two shot points (vee 
Suffield and Swift Current), and secondly, that these 
recording locations, when viewed as a whole, should give 


fairly uniform coverage along the refraction line. 
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16. 


Once we had actually arrived at the recording 
location specified for the days shooting, we ehecked the 
bocablonifer powerhiinestiatlicthere were anya ws usually drove 
a mile or two north or south of this location in order to 
find a location free of power lines. In general, it was 
poseibtel fo findasuchpasbocationyc! Hawi ngearrivedgarta 


suitable location we proceeded to lay out the geophone cables. 


Laying out of the cable was performed as Follower 
one cable reel was placed at the back of the truck, and 
while one crew member drove slowly along the road, the other 
puLlded the cable*into the ditch at the edge of the road, 
As each geophone takeout was reached, a seismometer was taken 
fromthe trucks, unclamped,tand dugointontghe ieround@arro a: 
bocation,.such,.thateits* top was gust: flushawith! thessugeace 
ofthe grounds It was found that this gave the maximum 
Signaleto noise ratio. It should besnoted that 40 feeteon 
slackscable wasApaidvout \atethe beginning. reThis!washdonesd. 
ensure that the spacing between geophones 6 and 7 was 
consistent with the spacing of the other geophones. This 
same procedure was followed for the other cable. Having Jaga 
out both cables, the truck was driven back to the centerven 


the spread and pre-shot checks were begun. 
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17. 


After the cables were plugged into the amplifier 
bank; the bine testoswitch was used. to test the. convinuiiy 
of the lines’ If thee resistances on’ each*feophone werer Touma 
to be in the range of 2700 to 3500 ohms, they were considered 
to be all wight... If the resistance.was infinite, cage 
usually implied an improper connection, and it was then 
checked by one of the crew members. A geophone whose 
resistance was too low was checked first for a shorted connec= 
tion, and if this. was. not.the.cause,. then, the, instrument wae 
replaced by the spare, Having completed. the. continua is 


check, the next problem was to set the chronometer. 


Using WWV as our time signal, the chronometer was 
started, then’ set°against this time signal by® means® of *a 
manual control. A short record was then made of the WWV 
signal and the chronometer marks on the recording camera. 
This record sufficed to calculate the error in the chronometer 
against WWV . Absolute time could now be read off the 
record, of the, chronometeromarks"toxan’ accuraty2ori about 7s 


milliseconds. 


The amplifiers were now switched on, and the traces 
of the seismometers observed. The ambient noise level 
tolerable was a deflection of up to 1/4" on the high gain 
traces. Normal operational settings of the external switches 
were 16 cps high cut filter, and either 0 or 6 db. attenuatuen. 


Line balancing was always used. With the modification of the 
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amplifiers, the normal high cut has been increased to 24 cps; 
and the normal attenuation is 18 db. or greater. This figure 
is the lowest yet used, and this was at an exceptionally 


quiet location, 


The system was then shut down until about 5 minutes 
before shot time. At this time, all instruments were 
started again. The WWV receiver signal was observed on the 
camera, the standard broadcast receiver was set to a pre- 
arranged broadcast. station and its flilvered, Pectitred 
signal observed on the camera, the chronometer (which remains 
running until the end of the operation) was rechecked against 
WWV, and the ambient noise level was rechecked. Any adjust- 
ments to the various Tilter and attenuation settings Wer. 
made at this time and were considered to be final. Av a@bGiue 
five seconds before the shot time, the camera was switched on, 
and a record of ground motion was made until 2 minutes after 
the pre-arranged shot time. The record was then developed 
washed, fixed, rewashed and then examined for arrivals if 
they were not visible on the galvanometer bank during actual 
recording. It was then labeled with appropriate informatie 
as to shot time, location, date, filter settings, attenuausar, 
etc. The instruments were now shut down, the galvanometer 
bank reclamped, the lids put back on the developing cans, 


the cables unplugged, and antennae were taken down. We then 
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drove to either end of the line, backed up the line, with one 
crew member sitting in the back of the truck reeling in the 
cable, stopping at each geophone location to pick Up TRE 
instrument, clamp it and replace it in the box in the truck. 
When both cables were picked up, they were fixed in the back 
of the truck, tied down (in order to keep them from moving 
too far in case of an accident}. All the other loose geaq 
was then packed into the van, and we were ready to proceed to 


phe next locatiom 


1.3 Spread Locations and Shot Information. 

Fig. 7 is a map of southern Alberta and Sso0uthwescuesa 
Saskatchewan. The dashes indicate the spread locations, 
and are drawn approximately to scale (spread lengths were 
2 miles). There are some 13 records recorded from the 
Suffield shot point, and 1/ records made from the . Wie 
Current shot point. Those’ spread Tocations west of surltveam 
were recorded from the Swift Current. shots. This exteueiee 
was made in order to obtain Tirst “arrival iInformatiouse 


the Pa PeTracuLons 


There is a gap in the middle of the line. THisiae 
is located in the Great Sand Hills of southern Saskatchewam 
We felt that since map coverage was poor in this area, and 
the surface material was sand, it would not afford a good re- 
cording location. Subsequent perusal of the area confirmed our 
SUSPICLONS. 

Table 1 gives information regarding shot datesyauume, 


distance from the shot point; and @size of the sca | 
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Spread No. Shot Date and Hour Distance (kms) Shot Size (lbs) 


1 O48. 62°. 10,00 Th OS 1000 
: Setee Eee tO, 00 89,28 1000 
S OO 63) 10,00 12767 1000 
4 Owe oe LO; peel as, 1000 
5 UOC og elspa is 135.82 1000 
6 > hell Olhaatad es paren le eg Ol 8, 155.33 1000 
ye Omi Orn 4 1, OO 165.18 1000 
8 eM"TO" 62 “10.00 169.94 1000 
9 O10 P62 15.00 Ueilge 1000 
10 Q-i1, “65° 10,00 Los, 52 1000 
shal Ot e238 15700 Ol. 1000 
T2 GOTT 63 10-00 205.76 1000 
a oO Ti G3. 15,06 213.89 1000 


Swift Current West Profile. 


Spread No. Shot Date and Hour Distance (kms) Shot Size (lbs)} 
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e.l Power Spectral Analyses. 

Some seven records from the Suffield East line 
were analyzed, while fourteen records from the Swift 
Current West line were analyzed. In this thesis, however, 
only the “sevén spectral analyses of the Suffield date are 
shown in fi@madS vand threesof, the Swiktt cCurvent tspecvral 
analysesiwint Bhen.O as Whe sereason wom sithe wherge’ numer eo. 
Suffield records is that they were considered to comprise 
the most uniiorm group vor udataz thecshdt size pemaraed 
Gonsvant atcl00G Lbs. and rthe rshot notes were toni wide 
once, Ucihneviwifit Cumrentiadata, “howevera anemnov tieea1, 
so uniform. o The shot size warred in most cases scant colten 


the same shot hole was used more than once. 


The computational techniques are described in 
Appendix 1, but it might be instructive to déscribe the 
method used in digitizing. All the records were digitized 
by nand. The section digitized in every case consisted of 
that part “or “the "record from about COO msécse, “befora sic 
First arrival to about 5 to 7 seconds after the firs; 
arrival. The digitizing interval was 10 msece. Careuwae 
taken to insure that the records were lain flat, Sco tha 


only linear trends were introduced into the data; such 
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trends being removed by the computer program before actual 
spectral analysis was begun. Only the first 5 to 7 secs. 
Was Cleitized because nearly all of the tnAformation “of 


Mivceres, Was “contained in this part or the’ record, 


e.ll Power Spectra of Some Suffield Records. 

Fig. 8 illustrates the Suffield power spectra. 
ieerevare .sOme important spoints ainuréegard ono thése sdata: 
The wmnesttamportant ils atherfact thatithere -appears sto! be 
little or no change in frequency of the main peak (5 cps) 
Wren stance giromuehedshot spoint re (This el senct*expected 
emce jaccording to the theory of wave propagation in layered 
media, there should be an attendant attenuation of high 
frequency signal with distance from the shot point (McDonal, 
Ameona et al (1958), Ricker (1944), Kovach, Lehner and 
Miller (1963) ). We note that any fine structure in the 
beokow 1s. veal, since the relative noise level 1s some 


2 or 3 magnitudes lower than the main peak. 


The raexisnicet important Beatune.of this sev of 
Gace 16 the fact that the main band of energy transmission 
is 5 cps, an unusually low frequency for such close-in 
refraction work. A more normal frequency is considered to 
be about 10 cps. We note in the diagram that the geologic 
environment of the shot holes is a heavy glacial boulder till, 


overlying arenaceous Cretaceous beds (Jones, 1963). The 
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boulder till is also riddled with gravel lenses. This 

type of shot hole environment will attenuate the high 
Preguencyol the oureoing explosive pulse, such that only 
Yery LOW irequencies remain. The resulting records have a 
very "loopy" appearance, and correlation from record to 
record 18 very difrrenivsinces there As=virtually no character 
to the record (Dr. E. R. Kanasewich, personal communication 
and Me. We GASmith, personal communication). This last 
(ecw Probew.Ly accounvs for the major uncertainties of the 
proposed .chusval Stmuctiwe, “simice the Suffield travel time 
Braph ie necgconsidercd. (oO, be as reliable as etiireereverse 


OVoOrL Le. 


iM Gteinteresting to note that there appears to be 
a peak at abouts WOwerps, which, on the average, seems to 
diminish with avstence. This is misleading, however, since 
hicespeak 16 gipu*feal in all cases, on some records being 
Sienel, walle tom others it is noise. This can easily. be 
Seen on the original records by noting the presence or absence 


of 10 cps enemgymbesore the first arrival. 


Although no data past 25 cps are shown, spectral 
ealeulations were carried out to 50 cps. On some of the 
power spectra calculated there appeared to be a strong peak 
at 40 cps: we consider this to be due to aliasing of a 60 


cps signal which is sometimes evident on these records. 
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APPROX. POWER AMPLITUDE 
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POWER SPECTRA OF SOME SUFFIELD 


SHOT RECORDS 
SHOT SIZE+=!I000 Ibs. 


SHOT TYPE= buried in 
large diometer holes of 50 ft. 


4x =205-76 kms 


GEOLOGICAL ENVIRONMENT: 


Heavy glacial boulder drift, overlying 
| arenaceous Cretaceous beds. Shot 
holes drilled into the side of a 
hill of low velocity moterial. 
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2.12 Power Spectra of Some Swift Current Records. 

Fig. 9 illustrates some Swift Current power spectra. 
It is immediately evident that the power spectra are much 
more complicated for this series of shots. Perhaps we 
should first note that the noise is some two or three 


magnitudes lower than the signal. 


ihe most important feature of thie set of data is 
the fact that the high frequency content has not been 
eopeuluated.as in the Surfield records. This ts due in main 
po,uwie Shou. nole environment. The holes. were drilled into 
a vive clay horizon of the Upper Cretaceous. This being 4 
fairly compacted medium, with a velocity of about 6,000 ft/sec, 
fue Trequency data is transmitted quite readily. We should 
note that shots number 17 and 23 were shot in holes used 
for a second time: there appears to be little attenuation 
of high frequency data, even though one would expect this 


to occur since the shot holes would be brecciated. 


If one examines the highest peak of the three 
records presented, there appears to be a frequency shift 
downward, as the distance from the shot point is increased. 
Whether this is due to attenuation of high frequency with 
distance, or to an increase in shot size, is not possible 
to determine with the present set of data. Further work, 


with a more uniform set of data,is required. 
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Another important point to note in these records 
is the multiplicity of peaks which are visible on the main 
peak. Interference of two seismic pulses may be 
the cause. of’ these 'ipeaks. or, as is more likely, they are a 
reflection of the fact thatuwaves,irayelling along different 
paths would tend to have different frequencies due to 
aqditferences jin, path lengstint. weVusbbaexexamination of the 
records tends to bear out this latter statement. The first 
break (considered in most cases to come from the sub-basement 
refractor), has a frequency of anywhere from 10 to 12 cps, 
while the head wave from the Moho appears to have a 


frequency of about 8 to 10 cps. See P. S. Veitsman (1957). 


Asjopposed polthe Suffileld—records,.theseismic 
records from Swift Current have very well defined character, 
especiallyijin the case of the M refraction. Correlation by 
character is quite possible in many cases, and the phase 
velocities of most] correlative arrivals have not been 
anterfered with 48 ja result of the low frequency character 
so evident in the Suffield records. Jf one accepts 
correlations made on such a basis, then we must accept the 


Swift Current travel-time graph as being well defined. 


We conclude, then, by stating that the frequency 
spectrum of the two sets of data are very much cdépendent on 


the shot hole location. In the case of the Suffield shots, 
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EXAMPLES OF POWER SPECTRA 
FROM SWIFT CURRENT SHOT RECORDS 


Shot size: # 27 =800 Ibs. 
#23 =600 Ibs. 
#17 =450/bs. 


Shot type: Buried in lOOft. holes 


#2 


7 


Ax = 285-06 kms 


# 23 


Ox 236-86 kms 


#17 


Ax =/13-84 kms 
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FREQUENCY (CPS) 
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GEOLOGICAL ENVIRONMENT: 


Shot holes penetrate oa blue 
clay horizon of the upper 
Cretaceous _ 
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the spectrum consisted of one peak at about 5 cps, with 
little or.no,change in frequency with distance. This was 
presumably a result of the shot holes being in a low velocity, 
unconsolidated medium. The Swift Current records, on the 
other hand, have a spectrum which has clearly defined peaks 
out to a Prequency of: about l2 cps. There appears to be 4 
Sa Gin the main peak to lower frequencies with distance 
from the shot point, but it is not known whether this is 

due to a change in shot size or is due to attenuation with 
Gistance. The fact that there is more high frequency data 

in the spectrum reflects the fact that the shot holes were 
drilled into a blue clay horizon, which is a fairly consolidated 


medium, having a velocity of about 6000 ft/sec. 


eeu AMpLICcUdes:. 

Although not directly related to power spectra, 
one should examine the records to see if there are any trends 
associated with amplitudes of the first arrivals. A plot 


ef Ghis“data“Vs"Go+bé Ffound’in*figs" 10. 


As is evident from this eraph, there is quite a 
large amount of scatter in the experimental results. The 
points 1 and 3 are the only pair of measurements corresponding 
to records from the same shot. There is a marked attenuation 
of amplitude with distance noted in these points, but since 


6ther pairs of records could not be plotted, it is impossible 
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to draw any conclusions from such fragmentary evidence. We 
see that standard inverse square and inverse cubic curves 
(constrained to pass through point 1} do not fit the data 
aural well, In fact... thera is very little pattern which 


emerges. from this curve. 


The reasons for this lack of pattern are perhaps 
two-told. it is felt that the shot, hole location and 
recording locations are responsible* for this Lack of 
coherency. As we pointed out earlier, this setvof data is 
cme mosc Unliorm set. This was e€videntyJéven tneéthe craters 
formed by the shots: their diameters varied only some 
sertee, in 30 freer. Although this, appears insignificant, if 
one considers the volume of such a crater (assuming spherical 
symmetry), variations of 3 feet in 30 feet correspond to 
variations of about 30% in the volume. There would be a 
concurrent loss in energy in such a shot; we feel that such 
large differences in the outgoing pulse energy would account 
to. a large extent forthe scatter in the amplitude measure- 
ments. Furthermore, variations could be expected in the 
geophone-to-ground coupling from location to location. 
Differences in the surface low velocity layer thickness 


would also account for differences in amplitudes. 


To sum up, the scatter in the data is thought to 
result from two sources. They are loss of energy as a result 
of brecciation of the’ shot hole and loss of energy due to 


variations in the low velocity layer thickness and geophone- 


to-ground coupling. 
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Chapter 3 


3.1 The Post Precambrian Surface Structure. 

From the travel time graphs of the experimental 
data we were able to deduce four separate refracting horizons. 
Well velocity logs, kindly supplied to us by the British 
American Oil Co., Ltd., enabled us to draw up a reasonable 
model Tor the néar suriace structure. These ‘lrogs defined 
Ouive clearly three distinct velocity aGistributions, a, low 
Velocity surface layer consisting of all post-Mississippian 
sediments pethe highrvelocity Mississippianselimestone; and 
another low velocity horizon, the Cambrian sediments. In 
order to simplify our model, it was decided that lumping 
in the low vetocitty Cambrian sediments with the surface 
material would not affect our model, and would leave us with 
two layers at the surface rather than three. Figs. 11 and 12 
{Wustrate Che post “rrecamorian sSuriace structure av Ssurfreld 
and Swift Current as derived from the well logs, while 
Figs. 13 and 14 show the modified models of these logs used 
ta our calculacrons. 

It is perhaps important to point out at this stage 
that these velocity logs served a further purpose. The 
experimental travel-time curves were quite complete in their 
information regarding all three refractor; below the Pre- 


Cambrian surface, but control of the Pe arrival from the 
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Low Velocity Layer 
NZ =. O74 kme. 
V 2 2. 51, ‘knis /eex.. 


2864" + .348 secs. 


Mississippian Horizon. 
DN Ereitcsinseny th Qo EN S4 
V = 5.12 kms/sec 


5615' ¢t 


Cambrian Sediments Sle secs. 


Az = .427 kms. V = 3.81 kms/sec. 


(O15* ¢ .624 secs. 


Precambrian V = 6.00 kms/sec. 
(Assumed) 


fell. §6O—6FdOosG Precambrian surface structure, as. defined 
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Low Veboci vy Layer 
Az = 1.138 kms 
V = 2.34 kms/sec 
83730" 
Mississippian Horizon. 
Aone fo ~KMS~. 


V = 5.32 kms/sec 


.486 secs. 


62071. -& .628 secs. 


Cambrian Sediments 
Az = .126 kms. “V = 2.42 kms/sec. 
an, O01 opt 
Precambrian V = 6.00 kms/sec. 
(Assumed) 


.680 secs. 


Fig. 12. Post Precambrian Surface Structure at Swift Current 
as defined by Continuous Sonic Velocity Logs. 
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Low Velocity Layer (+ Cambrian) 


V2 ,33ins/ see. 


Mississippian Horizon 
Az = .839 kms 
V = 5.12 kms/sec 
Precambrian V = 6.00 kms/sec 


(Assumed ) 


oes to... [Wo Layer Surface Structure at. curfield as 
derived ([fomeFileas bk: 


Low Velocity Layer (+ Cambrian) 
Az = 1.264 kms. 


V = 2.35 kms/sec. 


Mississippian Horizon 
Ne°= <755 (es 
V = 5.32 kms/sec. 


Precambrian V = 6.00 km/sec 
(Assumed) 


Fig. 14. Two Layer Surface Structure at Swift Current 
asnderived fromvfig. tle. 
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basement surface was not good in either direction. Assuming 
that the velocity of i was 6.00 kms/sec, theoretical 
intercept times”- were calculated’ for'this arrival at’ both 

Shot points using the well velocity data. This theoretical 
uime Was UCnéen used when fitting the Teast squares Tine to 

the basement data. The assumed velocity is borne out further 
by the fact that the velocities obtained from the lines 
fitted to the experimental data are very close to the 


assumed figure of 6.00 kms/sec. 


As is well known, much of Alberta and Saskatchewan 
is covered with a thin veneer of glacial material of one 
form or another. TALS was very evident along our refraction 
prerile > we Nave noL, however, made any correction for this 
surface low velocity layer. At most spread locations no 
mirormatcion as to velocities and thickness of this. layer 
aresoavaLllabplé, and it was further felt thav euch a correction 
would not enhance the data enough to warrant the expenditure 
ii GLMe . 

3.2 treparalion Of Travel Time Curves. 

See LLL 

The time of the shot break was marked on the spread 
Vocation records by one’ of’ séeveral methods.” As each truck 
was equipped with a chronometer, a WWV receiver and. & 


standard broadcast receiver, the signal from each of these 
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was recorded on the seismic record. Before a days operation 
the shooting truck and the recording trucks decided which 
broadcast station would be recorded. Once on location, the 
chronometer was set from the WWV time signal, and if the 
signal was strong enough, the chronometer was rated by 
Simultaneously recording the chronometer second marks with 
Ghose of the WWV time signal. iIf,; however; the signal°was 
Weak; then one -nad>to rely“on a@ different correlation’ device. 
Bown the recording trucks and the shot point truck recorded 
Give Tectitied, filtered, output of the. broadcast receiver sev 
Ho. wae pre-arranzed station. Then, if one did not have the 
WWY time signal to correlate by, one could correlate the 
shoo point and spread location records by using the 
prosicast® Station sliecneal. In most cases 10 was possible 
Go.merk the shot break on all the records to within 2 or 


Semi iii seconds . 


aeoees Distance .Caleulations 
The geographic coordinates of both the shot point 
and each recording location (usually one read the coordinates 
of either the first or the last geophone in each location) 
were readwoff a.large scale topographic.map...Anraccuracy p)of 
+ .05 minutes can be expected from such a map. These data 
were then put on punch cards and the actual distance calculation 
Was..carried out .on»an -1BM 1620, digital computer. The 
Department of Terrestrial Magnetism of the Carnegie Institute of 


Washington kindly lent the program to us. Hosmer'ts (1929) 


method was used in the program. 
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the phase velocities of each arrival in about the 
tives (ive seconds after the first arrival were calculated 
fer-each record. This was accomplished by plotting the time 
interval between the trough or peak of an event on each 
trace against the geophone spacing. One then fitted a 
straight line through the points and calculated the mean 


Melocity. An example of such a Plot i186 shown in figure 15. 


3.24 Travel Time Graphs 

Once all these operations had been carried out, one 
eould then begin to plot the travel-time curves. The type of 
plot used was the so-called reduced travel-time graph. 
This Variecy of lot was usede because one could expand the 


time Scave, thus enabling fone to make more accurate correlations. 


Hrom each arrival time on a particular record was 
subtracted a time equal to the distance divided by 6.5. 
Hence, arrivals falling on a line whose slope would normally 
be the inverse of 6.5 would now fall on a line whose slope 
was zero. Each separate arrival on each record was plotted 
Mreros cormecy position in. time and Space on the travel-cime 
Curve. hen using both the records themselves and the travel 
time curves, lines were drawn through arrivals which were 
thought to come from the same refracting horizon. Once we 
were gure that we had indeed made the correct correlation, 


ene Line was fitted to ite corresponding points in the least 
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squares sense. From these statistical analyses we were 

able to deduce the zero distance time intercept of each line, 
the apparent velocity of the refracting horizon, the standard 
deviation of the slope, and furthermore, we could establish 

the 90 % confidence limits of the slope and the Invercep. times, 
Both of these latter two quantities are, of course, minimum 
uncertainties, since there are a number of other errors for 
Wiich’We can-maké rio account. “These typesnof errors: biclude 
Ceviavlonsein thé-refracting horizonscornrom’ planarwaurfaces, 


hiecorrerationgor events; canlsotropy cimvythe horizon, evc. 


oe lewhous Of Correlaving Retracted Arrivals 
P. S. Veitsman (1957) outlines the various criteria 
used to correlate various wave-groups from record to record 

Oe arconvinuouslirefractilonsprorileszb) Lnvessencay these 

eri rerigcare fas -Tollowss 

1. As one progresses further and further from the shot point, 
the refracted head wave which is the first arrival will 
correspond to a refraction from successively deeper 
seismic °boundariés in 'théevearth's crust. 

2. The graph of the first arrival of each group of waves 
should be approximately rectilinear. Any departure from 
rectilinearity may be caused by relief on the refracting 
horizon or supraposed division boundaries. 

2, The apparent velocity of successive waves in the same 


group should be similar. 
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The reverse times along a particular refraction line 


Should be the same for any particular wave-group. 


The apparent velocities of one wave group should be 
Similar, whereas the apparent velocities of a different 


wave-group should be rather different. 


there should be a relative maximum in amplitude’ in’ 6ach 
Wave group, provided. that one is within the zone of inter- 


perences.' 1.6%, the zone where two dirfferent wave-croups 


Paperrere: due to .thei® similar arnriva. tamess 


There should be a determined regularity in the behaviour 
Gm the. amplitudesontene waveesroup Trom+location te 


Fogaptonr., aud Furr nermore, 1 rom prota teyn0 protise. 


There should be an inherent stability in the form of the 
proce. Aeliay vee, the waves oly one eroup rSshouwld tavesa 
Siuviar appearance [from record to, record, | We nave found 
Phateonts.ts particularly "(ruc o une Meng) retrace torn, 
especially in the region, close Vo yuhe cricvical disvance. 


Pom an example of Chis traitiseer ie. 16. 


Velctaman svaves that there should be a differences a tic 
frequency of the arrivals from different horizons. He 
cites as examples from his work that the predominant 
frequency of the P* wave-group was from 9-10 cps, whereas 


the P. wave-group had a dominant frequency range from 
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3.4 Discussion of the Travel-Time Graphs in the Light of 
these Criteria. 

Figs. 18 and 19 illustrate the Suffield and Swift 
Current Travel Time Graphs respectiveiy. sDurine theidiast 
Stages of correlation, the above criteria were used to-a 
large extent in making the final choices. A number of 
problems arose, however, one or two of which have not really 
beeneresolved satisfactorily; rather«we can only makewan 
attempt to explain part of the problem away. 

Criteria 1, 2, and 3 were satisfied on both travel- 
Gime curves. Criterion 4, however, is not satisfied. ‘The 
reverse times bevween the Sufiield and Swift Current shot 
points do not agree, the error being of the order of 400 
Mei Tseconas, oSoince there is so much data alone the Switt 
MUgtenu tine, we have decided to pus all our faith into this 
Picerprecaclon, and lump all €rrors into Une Ssurfileicd inver- 
Precacron. in Tact. we can show that this 18 nov an 


unreasonable assumption. 


We have seen that power spectral analyses of the 
Swift Current data indicate that the main bands of energy 
transmission are centred at about 5 and 10 cps, while 
Siigier analyses of the Suffleld data indicate that the 
dominant frequency is about 5 cps. Let us assume that the 
velocity of propagation in the earth is not frequency depen- 
dent; hence the reverse times along a path should not depend 
Griclcally on the frequency. This. tact is. provably true since 


the frequencies of interest are low. 
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Having established that there is a fundamental 
difference in frequencies along the reverse paths, we can 
account for at least part of the time delay. First, the 
shot noles at Suffield were drilled into a rather heavy 
boulder till, whose velocity of propagation is low. In 
addition to this, Jones (1963) points that there are a 
number of gravel lenses to be found in the general area of 
the shot point. Gravel beds are known to break up a wave- 
front rather badly, and have very low propagation velocities. 
such an environment would certainly delay the wave front, 
although the magnitude of this delay time is indeterminate. 

A further. consideration is that the resulting records from 
such a shot point have a very "loopy" appearance, there 

Demi sitcle or no trace of character from record to. record, 
Under such conditions correlation of events is very difficult, 
and the validity of interpretations made on such events 

is questionable. A further instrumental time delay is caused 
by the low frequency content of the signal. We should point 
Car. that the high = cut filter settings for most of the 
Surrielo shots was 12 cps, while that of the Swift Current 
shots was 16 cps. There is a time delay of a few milliseconds 
due to this filter setting difference. Last, and perhaps 

most important, is the fact that. the short polncervand: Ggeveclor 


locations were not exactly reversed. Thus the reverse records 
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could not be tied exactly. Another factor which complicated 
the tying of events lay in the difference of the frequency 
content of the two records. There was no similarity between 
any of the events on the two records which should most 


nearly reverse time tie. 


Criteria 5 and 6 were found to be true in general, 
but phase velocities of presumably similar events on the 
recerdss«eirom the Suffield line.did not,always.match..,This 
is considered to be a result of the "loopiness" of the 
Pecoccs.and.aringing,in the filter system, with a consequent 


rune Goeether of the various arrivals. 


Criteria 7, 8, and 9 have been discussed previously. 
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Chapter 4 


Here we shall present the final results of the 
calculations carried out to arrive at the Droposed crustal 
model. The arguments both for and against will be presented, 
and the reasons for adopting the final model even in shslad 2, lila 
of Une argument: against it. This model will be compared 
to that proposed by Weaver for the crustal structure west 
of Suffield. Finally, we shall present the results of the 
500 ton trial at Suffield and give a tentative interpre- 


tation of the results. 


4.1 The Layered Earth Model for Saskatchewan. 

Using both continuous sonic velocity logs and the 
travel-Cime graphs derived from our experimental results, 
we have decided on a six layer earth model. These layers 
have been called the low velocity layer, the Mississippian 
horizon, the Precambrian basement, the sub-basement 
retractor, the Conrad discontinuity, and the Mohorovitié 
@i@scontinuity. The final proposed model®*for-tnhe crustal 
structure in Southern.Saskatchewan may be seen in fig. 1/7. 
Below is presented a set of tables which give all information 
regarding the various parameters used in the calculations. 


Where the errors are presented, these represent the 90% 


confidence limits. 
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DEPTH IN KM. 


PROPOSED CRUSTAL STRUCTURE 
IN SOUTHERN SASKATCHEWAN 


Suffield Swift Current 


LOW VELOCITY LAYER V2#=259 km/sec. 


MISSISSIPPIAN LIMESTONE Vz" 5°22 km/sec 


PRECAMBRIAN BASEMENT Vaz 6-00 km/sec 


6 


= 19'E. 


20 
SUB— BASEMENT REFRACTOR V=6:51 km/sec. 


25 
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CONRAD DISCONTINUITY V2 7 253 km /sec. 
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MOHOROVICIC DISCONTINUITY V,8-OI km/sec. 
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47, 
West Shot Point (Suffield) 


Horizon Velosi Gy Thickness Intercept Time 
(kms/sec) (kms ) (secs) 
Tow Velocity Layer 255° Lt. 30 
Mississippian Horizon pau . 84 
Precambrian Basement 5.977 wee oto 
Sub-basement Refractor Ox47 ty. 08 26.15 2.031 = .166 
Conrad Discontinuity Teel. £4.08 19.15 6.068 + ,263 
Mohorovicice Discontinuity 8.06 + .11 9.419 = 416 
East Shot Point (Swift Current) 
HOrlZon We lool Gy Thickness Intercept Time 
(kms/sec) (kms ) (secs) 
Low Velocity Layer 2.35" 1.26 
Mississippian Horizon oe 76 
Precambrian Basement 6.00 (assumed) 8.87 im hoe 
Sub-basement refractor 6.55 ¢ .10' § Mh br “es. 36@s 725 
Conrad Discontinuity (eo Re 6.86 6.093 + .060 
Mohorovicic Diescontinuty 7.9972 ..05 O.g22 = 47 


The travel-time curves from which these data were derived 
may be seen in figs. 18 and 19. The other two parameters of 


interest in the final model are the true velocities of the 


* These velocities are calculated from the continuous sonic 
Velocity logs 


+ These velocities, and the velocities below this refractor 
represent apparent velocities. 


' This is a theoretical intercept time calculated from 
well log data. 
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layers and the dip of the layers. They are as follows: 


nord zon 1 i: frue VeLloc icy 
(degrees) (kms/sec) 
Low Velocity Layer O 2.59 (average) 
Mississippian Limestone © 5. 2eetaverace ) 
Precambrian Basement O 6.00 (assumed) 
Sub-basement Refractor LOtar 6.51 
Conrad Discontinuity 16' W 7e28 
Mohorovicic Discontinuity 2.4 Ww 6ron 


4,2 The Travel-Time Graphs. 


Thestravel-timescurves showneinnfizgsspis"ands19 
represent yunesdatarfrom ther Suffieldsand Swift. Current shot 
peatea respectively.) SomenlO-records weme tsed torpcompile 
the Suffield travel-time graph, while 15 records were used 
Poncompr lest ne Swikt Cunrent secraphs aThéese vecurvescwere 
duectissecdecariter, (but wenmustunowediscussothemainoeterms of 
whenproposedacrustal «structure. «Inrbothecasés we note lthat 
MiehiaLencepL tameschave kbeenlconrrectedovovthe first opreak, 
phe correlations having -beenamadenin *theaefirst troughs 
Arrival amplitudes have been graded for the Swift Current 
G@eabae woulnotefdrttheisuffireldedatasypstncevtneclatterseane 


nov considered tio be asrreliable: 
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feelin curtiel d.Shot sPoint. 

Lhe only .worigon, for which .cood first arrival 
control exists on the Suffield curves is the sub-basement 
Po itaclom., slbere bs little doubtvaste,the validity of 
this correlation. The basement refractor has been con- 
strained to fit a theoretical intercept time calculated from 
Beet loeces thesnumber of points used to fit this line are 
Mmoninal.  Gorrelabion..of,the.Conrad,and Moho .events<is not 
Sent dered to be .wery reliable inthe, light of the,low 
Preovcncyecharacter,of the» records...This isualso.evident in 
Thesstavistical errors calculated for the, slope and inter- 
Sete vines parameters. —The.error.bn,the.invercept times.is 
eppeciallyehigh for both these.events; infact, 1Lt.is,high 
Shouchatboseaccount: for. the non-reverse,time.tlie.»,spread 
velocities are not everywhere consistent with the calculated 
eoparens velocities. This ,is,especiaily evident.on._the 
sub-basement refractor line and the Moho correlations. 
fuchaditsacrepancies:.tend,to.further.reduce.tne.reliability 
SOfnune. data... We should keep. in_mind, .however,. the, low 
frequency character of the data and the sharp band-pass of 
Phoeuy?iiieraausea for most.of the,.recordings. ,Both,of .chese 
facts would tend to give erroneous spread velocities since 
the events are not clearly separated from one another, with 


consequent interference taking place. Although we. cannot 
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Glaimsinverference,forathe,first.arrivaly weade-note. that 

a high spread velocity is sometimes consistent for all 
arrivals throughout the record, while on other records this 
LS bOLpthe.case.».No.good,explanatilon+is»available forsthis 


phenomenon. 


4,22 Swift Current Shot Point. 

The Swift Current travel time graph, on the other 
hand, is much more reliable in many respects. We see from 
BugesbOn.vhnav,.the .first arrival .control fer two out efathe 
four refractors is very good (the sub-basement and Moho), 
while scontrol.tor,the basement .refraction.is.much,) better 
than that at Suffield. Once again the basement data have 
beenpconstrained,.to;-pass through.a.theoretical intercept 
LiMentere iS evident onathe travelL«~time,.curve, that imore 
tnan.one correlation is.=possible for the: basement-event.. 

We have chosen the final correlation for the following 
peacsons... «kirst.,.,bhe apparent.veLlocities,af the-«arrivals 

used for this correlation were quite consistent with the 
assumed basement velocity.  Secondly,,.the .close-in,.refraction 
work (records 13 and 15) was done subsequent to the original 
pick, and the data from these two records indicated that 

Ghe f£inal.pick .was.more reasonable.» .lhere ais .~ertaLnly «ao 
doubt about the sub-basement and Moho arrival correlations, 


since in both cases correlation by character was possible; 
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spread velocities agree very closely with the apparent 
velocity calculated; and first arrival control was VEry@o00 . 
In order to illustrate the definitive nature of correlation 
by character we can consider the M arrivals picked on records 
18 and 19. These arrivals were very strong, and had a 
character which is considered to be diagnostic of this 
arrival (see fig. 16). Correlation of the events on 

records 21 and 23 with the former events gavé a cood' Tine, 
whose velocity was close to that defined for the Moho. 
subsequent refraction work, (records .22,.\24, 25, 26.and 2%) 
Deve food tirst arrival control on, the M.refraction,. and va 
line drawn through these events tied almost exactly into 
Bieeti.ne Orawn on the basis of second arrival inrormation. 
Furthermore, the statistical analyses tend to bear out .the 
Correlations we ave made...The Conrad event 2S. Teli to spe 
G@iie well defined. Statistical calculations made, sub- 
stantiate the rectilinear nature of the correlations, and 

Gne time intercept error is especially low... There is 

always a certain amount of doubt in using second arrival 
frlormavion, but the doubt in this case is considered minimal. 
Having discussed the travel time curves as separate entities 
aavorveed from the crustal structure, lev. us now consider the 


two curves in conjunction with the proposed model seen in 
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feoeee vos and Cons of the Proposed Crustal Structure. 
Hes) WCong 
There are two serious objections to the proposed 
erustal structure. These are the non-reversing time tie, and 
the fact that the apparent velocities for the M refraction 
Beesnou consistent: with the dip shown,in, fie. 17.. The, dip 
shown is not calculated using the usual equations, but is 


Simply the dip as defined by the depths under each shot point. 


The non-reversing time was discussed previously. 
Aicnoughn we could not) account. frorcitefully,#wervecould ‘ar least 
Paola parvocof the*discrepancty.s .dWerturther novteythavs the 
errors’ in the intercept times’ for thetvarious’ arrivals: at 
Suffield data are almost high enough to cover the 400 milli- 
eecondsfecrror in the time tie. \eWeeteeli that thistis not a 
very serious objection in light of the poor data along the 


put Lela line. 


Mie fact that the Apparent Velocities are”*very nearly 
equal is a more serious objection, this condition implying 
little or no dip. The velocity under Swift Current (7.99 kms/ 
Bec) 18 consistent, while the velocity under Suffield 
(8.06 kms/sec) is too low to account for the dip indicated by 
the intercept times. The errors nearly cover this uncertainty, 
although this does not afford the best explanation. An 


alternative explanation is 
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available. If we assume that the dip becomes shallower under 
pulrlela, then we can easily account for the close apparent 
Welocities. Weaver's model does, in fact, suggest that the 
Gro uraer ourfleld 1s much lesa than that which we have 
Pieeaged Ii tie. 1, If this 1s °itideed the case, then 
We Meeca NOt worry about the discrepancy in the velocities, 
since this latter modification would justify the observed 
Vebules . 

ihe fined op jecvion to the model lies in the lack 
of coherency between the model proposed by Weaver (1962) 
rox che structure between Suffield and Vulcan, and that 
proposed by us. In essence, the two models are similar in 
that they both are six layer models, and the crustal 
thicknesses calculated are quite compatible. We do not 
Pestudate a low velocity layer as Weaver has done since we 
heave wo cate to indicate that such a layer exists in 
Saskatchewan. Although Weaver substantiates his argument 
for such a layer very well, we see no reason to include it 
1n-our model. The reasons are that we have no discrepancy 
between. the theoretical critical distance for the M arrival 
and the observed critical distance, and furthermore, the 
Bouguer gravity anomaly suggests that the low velocity layer 
pinches out just east of Suffield. There is some suggestion 
of this layer observed on the Swift Current uravel-time curve. 
The last three Moho refractions appear to come in just a 


little late suggesting a delay for these locations caused by 


a low velocity horizon (Fig. 19). 
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4.32 Pros. 

There are several facts which we feel substantiate 
our model, and most important, override the arguments 
against it. Let us first consider the Suffield data. 
Although we have already pointed out the many discrepancies, 
the intercept times calculated for the Moho and Conrad 
events agree closely with those calculated by Weaver (1962). 
Secondly, from the extension of the Suffield west data into 
the mountains this last summer, the intercept time cal- 
culated for M at Suffield was about 9.25 seconds; again 
agreeing very closely with what we have calculated. This 
is an indication that the picks made are not completely 


miscorrelated events. 


Certainly the exceltent’ first’ arrival data shown 
in the Swift Current curves substantiates the structure under 
mace sne. pornt: “The excellence’ of the* Statistics further 


bears witness to the reliability of this profile. 


Even though this model represents a straight- 
forward interpretation of the data, it fits extremely well 
with the observed data. The Moho critical distance cal- 
culated for it is very close to that observed on the records. 
We have calculated a critical distance of 132 kms. for the 
model. The observed critical distance (as judged by the 
location of the first strong arrivals from M) on the graphs 
ig about 125 kms. for the Suffield profile, and about 135 kms. 


for the Swift Current profile. This is a remarkably good fit, 


and helps to strengthen our interpretation. 
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4,4 Comparisons with other data 


1) McCamy and Meyer (1964) have shown a number 
of crustal profiles completed in Montana, Jue tC. SCUGh of our 
own line (see fig. 20). The most important of the profiles 
is in Eastern Montana.. .The/model they obtain for the crust 
Pe citiilar in allwrespédts/to that \pestwlated by us. 
Although the near surface structure is one layer of 2.8 
kms/sec, there are four other layers which correspond to 
our Precambrian basement, the sub-basement refractor, the 
Conrad discontinuity and the Mohorovitic discontinuity. 

The velocities for each of ‘these horizons are 6.1, 6.4, 

me and 6.3 kms/sec respectively "/-Thée dips \snown fin their 
diagram correspond in general with those found by us. 

iiao Such close correlatiom should exist between Lwo. profiles 
done completely independently is a good, indication that the 


model for Saskatchewan is generally correct. 


2) Weaver (1962) discusses the work done by 
Richards and Walker (1959), and has replotted the data used 
by them; arriving at a new interpretation using rather 
different results. His final model suggests that the sub- 
crustal velocity is about 8.2 kms/sec and the Conrad velocity 
vee avout («2 kms/sec, the crustal thickness being about 43 kms 
and the depth to the Conrad being about 26 kms. The interpre- 
tation used by Richards and Walker suggest the same velocities 
but the depth to the Conrad was 29 kms (as opposed to 26 kms ) 


and the crustal thickness is about 43 kms (the same as Weaver 
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obtained). The models obtained by Weaver and by Richards 
and: Walker are not at. all at ‘odds with ‘each outer. ite act, 
their relative closeness would suggest that although the 
authors are not in agreement with each other, their 
separate interpretations are not invalid. The present 
author has not recalculated the work done by Richards and 
Walker, but rather wishes only to point out that the depths 
to the various horizons postulated by them are shallower 
than those given by both Weaver and ourselves. A very short 
profile was shot (only 132 kms.), and scanty data were 
obtained from the records. The interpretation was based on 
second arrival information for the Conrad and Mohoroviéié 
Greconvinuities, which at the best of times 18 not nearly 
so g@00d as using first arrival information. Furthermore, 
1ncverpretation of seismic data is often subject to personal 
biases on the part of the seismologist, this being a bias 
Werehn 18 hard to account for. Reinterpretation of data is 
hacurally £oing to be subject to a different. bias and the 
resulting interpretation cannot be expected to be exactly 
the same as the original. Perhaps a more accurate picture 
would be a composite of the two. We do not feél that the 
Richards and Walker profile is completely in error, rather 
we feel that the profile is too short to allow a very 


detailed description of the crustal structure in the area. 


i? | ve 

abrado Lt vd bos tevesW ¥0 ‘sabe Bet en | es = 

Jost of .netite nope Abe Shun te ifs oe fon ' 

si? ravoditie ted? dessada Bittow esoae4o ts ovtse a ben terit 
vigod .1edto tose Athy dnseeetas 

inegetq ol? .bEleval Somers venoisesonayezat oo 

brs ubtatoth yd 506 dow ont retatitiaoes oul | 


jan? duo galeq od yl 


ce 
- 
5 


towollada ous madd yd bSvsLaseog ands ld Prone 
jvone yrev A .asvfeote bose teveeW dod eet at 
aitow stsb yvinsoe bate . ts ani Sef vieo) tors « : 

tn 

no) Bsese sew ooldavetqesdnl sar Lebirooet etd pres rf 


Si5lvouccdtoW. bas. hetaed said a07 pobienotat Levi 
'yson gon at semtt Yo teed ese dB co 00 rt 


- oy — 7 a 


stomata .aotyantotat fev ies sane ile 8B boog 
lercarsg o¢ Jos@loua ostto @h sdab o feted to 10. ete 


y 


aid B acted ated? .tebgelomagee sad To. rag oF 
bi sisb te solssisiqvednisA “tot Javoo98 om 
sciJj bres este Jnerettib » ot sootdus od of 3 
yidosxa os o¢ betosqxe sd Fonts noted [rou 
TU90 5q evan unos etom & aqantet eine ro: az | 
ait tadt Lest con Ob sw. — Pr to. Mais no 2° 
rast sorte ot LStoLqaeD at - ji _~ W bes | 

yrtete 6 wot 16 ne AGS marie 
esta ond AE baubouon aithy ; 


60. 


Although the crust proposed by them is thinner than the 
feelns Suructure, it is in line with data recently obtained 
pyrourselves in the profile shot from Suffield west into the 
Rocky Mountains. These data suggest that the crust thins 
under the foothills, the velocity being about 8.2 kms/sec. 

A remarkably close correlation thus exists between the data 
obtained by Richards and Walker and our own data, even 
though the Rocky Mountain data cannot be subject to a 


definitive interpretation (see 4.5). 


3) Recently Hall and Brisbin (1961) completed a 
Bieie Cnoed Telrracviol Droliie near thew Flin Flon area if 
Meantvoba. ~The principal results obtained from this survey 


ehiget 


a. A three layer earth was postulated: 


16.5 £°1:1*kmes 


Il 


depth to the Conrad discontinuity 


beour to une Mono GUisconvanulty” “= 55,2 2 le" hme 


These values are markedly different from those obtained by 

us. We should, however, point out that the seismic-refraction 
survey was carried out along or near the boundary between 

the Churchill and Superior geologic provinces. Since these 
two provinces show such different ages, it is not unreasonable 
to assume that there might be quite different cOonaLcLons 
existing for the deep crustal layers. The Moho dip calculated 


for the Saskatchewan data is such that if it is extrapolated 
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into Manitoba, the crustal thickness would compare favorably 
with that described above. The Conrad in Saskatchewan shows 
PetierKablyalitole di pyacand ds mot «taal compatible with that 
described by Hall and Brisbin. The fact that the Manitoba 
profile was shot near the boundary of the two provinces is 


fell influence this norizon..to..some extent. 


ime Pheu velocities are ‘as follows: 


L2G. 
I 


6.15 + 0.05 kms/sec 


Pp! =! 7,10) £0.05) kms/see 


a0 
i 


8.17 + 0.06 kms/sec 


These velocities are well within the range of the accepted 
parameters, and agree quite closely with those obtained by 


ourselves. 


4) Evidence for intrabasement discontinuities is 
found in many areas. Robertson (1963) has found reflections 
from an intrabasement layer in Southern Alberta, in an area 
south of Weaverts line. Work by Junger (1951 ) it) Base Horn 
County, Montana and by Widess and Taylor (1959) near the 
Wichita Mountains of Oklahoma has indicated that intrabase- 
ment layers exist in other areas. A similar discontinuity in 
our own area might be the source of our 6.5 km/sec arrival. 
imilicht ofthe complicated nature’ ‘of the’ Precambrian base- 
ment, the existence of such layers should not be too surprising. 
They might represent the roof of a batholith or some other 


large intrusion, or they might represent the boundary between 


different geologic provinces. 
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We do not have enough data to substantiate any of these 


hypotheses. 


AS a consequence of the foregoing remarks, we feel 
that the model postulated for the earth's crust in southern 
Saskatchewan is a reasonable model. The aereement wt the 
theoretical model with the experimental data is remarkably 
good. The fact that an independent refraction Prot le .shiot 
in Montana, just south of our line, agrees in Benes l orn 
with our own is further justification for our model. The 
statistical analyses of our data are certainly good for 
the Swift Current line, and the layered structure at this 
shot point is believed to be correct. The agreement of 
crustal thickness at Suffield between our own model and 
that of Weaver for the western plains is also very good. 
Alinasthese,tacts increase our .confidence in the preposed 


sbrucrure. 


4.5 Suffield 500 Ton Trial and the Rocky Mountain Data 
In July of 1964, the Defence Research Board of 
Canada, at the Suffield Experimental Station, detonated 
500 tons of TNT on the surface of the ground. Ten seismic 
records were obtained from this shot, ranging in distance 
from 200 kms .to..500 kms. from the shot point... Bight of 
these records were provided through the courtesy of the 
eae coce ss. Cmistal .Sbudi 48 Group seismic crews. one was 


supplied by Dr. Hugh White of the Dominion Astrophysical 
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Observatory, Victoria, British Columbia, and one record was 
provided by ourselves. These results were compiled and 


plotted on a reduced travel-time curve (Pie, sey. 


This diagram illustrates first the surface tTopo- 
graphy, and below, the reduced travel-time data. The reduc- 
tion factor is X/8.25. This reduction factor was chosen 
Since Weaver obtained a velocity of 8.25 kms/sec for the 
Moho under the plains. The squares and the triangles 
represent the reversed travel-time data used by Weaver, the 
anmvercept time of this data being about 9-2 secs. This 
time corresponds to a depth of about 46 kms. The circles 
represent the new extension of the lines, also reduced at 
ope kms/sec. As one readily observes, there_is a 
aerinite structural implication in the reduced times, We 
musc, however, point ouvetvhat "these data are in the raw 
Svaue, Salice no correction has becr™appiiecu.o accoun, for 
eLevarzonwor for the high speed overpurdew which existe in 
the mountains. The effeets of these factors are discussed 
below. Since no reverse profile was shot, a definitive 
interpretation is not possible; as a consequence, we shall 


present several alternative interpretations. 


All the points have been displaced to the point of 
their approximate emergence from the Moho; thus any structure 


in the time domain is in its correct location with respect to 
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the surface topography. The first interpretation is the 
following. Let us assume that the depth of the Moho remains 
constant under the mountains (this would be about 46 kms as 
derived by Weaver). The changes in slope indicated on the 
Grave l-time eraph then imply "that? there Is @ VELOCITY 
gradient from east to west starting under the foothills. 

ine change ity slope’! is equivalent® to a velocity OF 8.5 kms/sec 
wider the foothills (see Tig. Cl) breaking to a velocity of 


Ont kms/sec under the Rocky Mountains proper. 


The next inverpretation is one in which 1t is 
assumed that the time YVanuetvions are due to structural 
changes under the mountains, and not to changes in the 
velocity),of the Moho. In this, case,. the. Moho appears to 
become shallower under the foothills, the order of magnitude 
being about 5 kms, and then thickening again slowly to the 


west under the Rocky Mountains. 


The Tinal, “and most probable inverpretvatvion,”* 15 4 
combination of the above two. In this case we would postu- 
late a velocity eradient. from 8.25 kms/sec im the, plains: to 
about 8.00 kms/sec under the Rocky Mountains. In the light 
of this velocity variation, there would be a thinning of 
the crust under the foothills, and then a roughly constant 
thickness under the mountains, the depth to the Moho being 
about 35-40 kms. This latter interpretation is in good 


agreement with that postulated for the Southern Rocky 


- a 
/~: i 
aa i? re _ dg 
ee a ue _ ; 
AOorn ; die a ah if 
sig af nolsssorgisdnt te art" uns qa oo 
ofteameat oom arit te doqeb stig ey omueeg. wat £ 
- oe 7 
as em oc? guods od biwow a tday! entavsuon orld ¢ 
“iz oo heteolbal eqoie Ak eaanaito. ad? ine 


pei 


vJloolev s ef ered feng yk: fond gory 5 por x oy te) 
eillivoo!) sad} sebim autoneda seaw or vase “hom on 
oeg\ann 2.8 Yo ydhooLav «6 o¢ jnstey Lupe aL ogole wae see Bigs 
to ytisofev « oF gntWweetd (£5 .git wee ektnavout * mie Th 


Ja 
regeTd amtstnwem ywilodR' ely soba ose 
= 


dotdw af ano el notistenqretnt oxet odin 
latuJouTde of awb ote enottedeey | emths oud re TL. 
snz3 oi sesosdo of Jon bre entetauon t: } is hats ast 18t 


2 
) 


w 
tt. 
— 


ot atasqde ofcM aid yses0 add’ of 


2nisdnvom ylooh uh ws 2 


s el ,coljstergzetial sidsdorg sOar. bas | ont? 9 ar 


-%e0q bluow ow sees eft? nt ont vod oa ute z ori. 
| : ok 


od atiieig eft ot cea\aqt 29.8 goats 


tdatil edi ol § .é@hiednvom ——- z) 


itetacoo eeplausies B ast biciae.) | 
afitad ovtolt., oid og ei sary Sit “a 
boow ae’ et pipalibians 
eit mei a 

A ; = J mS 


in 


ve 
2G ee 


Mountains (Mr. L. C. Pakiser, personal communication), 
the velocity gradient being in agreement with results 


published by Stuart, Rolleret ‘al (1964). 


We prefer this datter interpretation, but a final 
answer is dependent on the measurement of the true velocity 
of tne Mohorovicic discontinuity by means Of a properly 
reversed profile in the mountains. Below we present a few 
speculations regarding isostatic compensation and a dis- 
cussion of the effect of the overburden on the observed 


Mivercept times. 


4.51 The Effect of Overburden. 

Calculations made to correct for elevation and 
overburden in the Rocky Mountains indicated that the 
correction to the data is negligible. At most, the arrivals 
are some 0.10 seconds too early with respect to the arrivals 
in the plains. The following is a list of the residual 
time to be added to the intercept time to correct it to 


the approximate level of the plains: 


From west to east the time residual (in seconds) 
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By observing the curve shown in fig. 21 we see that these 
values do not change the shape of the curve, only lower 

ius west end by a factor of about .10 seconds in time. This 
only implies that the depth to the Moho is changed by about 
1 km at most (i.e., the Moho becomes 1 km deeper). Our 


final conclusions regarding the structure still stand. 


4,52 Isostatic Compensation. 
Although the data are unreversed, and not really 

COrrecvea TOY’ any Sory of crustal varlations, they do suggest 
that the crust is thinning under the mountains. We felt that 
ie wou.g ve instructive to look ‘at the two theories of isostasy, 
and do some calculations upon the data. The two most widely 
quoted theories of isostasy are the Pratt-Hayford Isostatic 
System and the Airy-Heiskanen Isostatic System. The 
assumptions of the Pratt-Hayford System are as follows: 

1.  idostatic Compensation is. uniform; 1.e., Une deneivy 
underneath mountains is uniformly less than under flat land. 

2. The compensating layer is located directly under- 
neath the mountains and reaches to the depth of compensa- 
tion D, where equilibrium prevails. 

3. The density p' of the compensation, corresponding 
fo a topographic elevation h, and the density p of the topo- 


graphy satisfy the equation 
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4. The depth of compensation is everywhere equal when 
measured from the physical surface of the earth and not 


from sea level. 


The assumptions of the Airy-Heiskanen system can be summarized 
as follows: 

1. Isostatic compensation is complete. 

2. The compensating layers are directly under the 
topography. 

3. The density of the earth's crust is everywhere 2.6/7. 

4, The density of the underlayer is also constant 
everywhere, having a density which is 0,6 higher than the 


density of the crust. 


These assumptions imply a root of crustal material 


under the mountains. 


From our data it is apparent that the Pratt- 
Hayford theory of isostasy is more applicable. Considering 
this as our isostatic model, we calculated the depth of 
compensation for the Rocky Mountains using several 
different models to demonstrate how this depth varies with 
varying parameters. We have, of course, assumed that 
isostatic equilibrium prevails in the mountains, but this 
4s approximately true (see Heiskanen and Vening Meinesz, 


1958, Garland and Tanner, 1957). 
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Below we present some of the various models we 
have used to calculate the depth of compensation. The 
following parameters were varied in the calculations: 
the density of the crustal material, the density of the mantle 
Moaverlal, and the average height of the Rocky Mountains. The 
densities were deduced from curves published by Nafe and 
Pie (1Ol?). ” Fig. “oo ta the schematre model used to 


calculate the level of compensation. 


Rocky Mountains 


Plains 


Mantle 


Mantle 


Level of Compensation 


PL ee 


The following parameters were kept constant throughout the 
calculations. The thickness of post Precambrian materials 
under the plains was 2 kms, and the density was 2.67. The 
thickness of the crust under the plains was 45 kms, and the 


thickness of the crust under the mountains was 40 kms (not 


including the sedimentary section). 
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Reference to Table 2 gives the various models 
calculated. We see that the average sort of model has a 
depth of compensation of about 116 kms to 118 kms. We 
further note that this level of compensation is rather 
insensitive to variations in crustal velocities and 
variations in the average height of the mountains. They are, 
however, very sensitive to changes in the mantle velocity. 
The change under the mountains of 8.00 to 7.80 kms/sec 
produces a change in the level of compensation of some 
35 kms. The numbers we have derived for the level of 
compensation of the Rocky Mountains bear a startling simi- 
larity to the level of isostatic compensation that Hayford 
and Bowie calculated for the United States, namely 113.7 kms. 
Whether our numbers are coincidental or whether they are 


significant remains to be seen. 


4,6 Conclusions 

From all the calculations made, we have derived 
a six layer earth model for southern Saskatchewan. These 
layers we have called the Low Velocity Layer, the Mississ- 
ippian Layer, the Precambrian Basement, the Sub-basement 
refractor, the Conrad discontinuity, and the Mohorovitie 


discontinuity. The velocities calculated for these layers are 


V, = 2.59 kms/sec Vy = 6.51 kms/sec 
Vp = 5,22 kms/sec Ves = 7.23 kms/sec 
V 6.00 kms/sec (assumed) Ve = 8.01 kms/sec 
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The depths under the shot points of each of the layers are 


as follows: 


Under West Shot Point Under East Shot Point 

Mississippian = 1.30 kms Mississippian = 1.26 kms 
Precambrian = 2.14 kms Precambrian = 2.02 kms 
Sub-basement = 9.91 kms Sub-basement =10.89 kms 
Conrad =36.06 kms Conrad =35.50 kms 
Moho =51.21 kms Moho =42,36 kms 


The dips for each of the layers may be observed in fig. 17. 


Although this model does not agree with that 
postulated by Weaver for the crust just west of Suffield, 
there are many features of our model which suggest that the 
true structure to the west is not the same as that which we 
have found. First, the Bouger anomaly which Weaver uses 
as-part of his proof dies out just east of Suffield. 
Secondly, we have calculated the critical distance for the 
Moho refraction and find that it agrees very closely with 
that observed on the records (about 130 kms). Thirdly, 
the crustal structure which we have postulated agrees very 
well with the crustal structure in Montana, along a line 
parallel to our own, but south of us (see McCamy and Meyer 
as well as fig. 20). Although the Conrad does not dip very 


much, the Moho dips some 2 degrees to the west, and the 
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implication is that eventually the Conrad will pinch out to 
the east. This agrees with present data in eastern Canada 
although the Lake Superior Experiment may reveal this layer 
(Mr. Malcolm Bancroft, personal communication). The sub- 
basement refractor which we have found is not unreasonable 
since other people have reported seeing reflections from 


intra-basement layers. On the whole, we feel that the model 


which we have postulated is a very realistic model. 


The extension of the Suffield data to the west 
Nas produced some interesting results. ‘Until either ‘converted 
Wiapes are analyzed, or a reverse profile’ Te snov, 2 
definitive interpretation is not possible. Some alternatives 
have been presented. The most reasonable interpretation 
seems to indicate that there is a west-east velocity gradient 
of about 8.25 kms/sec under the plains to about 8.00 kms/sec 
underasthe Lodthidtwe.n( in addbtion, dtinippearsmihawanhe-cruist 
thins from 46 kms under the plains to about 40 kms under the 
foothills, and remains constant at this latter value under 
the Rocky Mountains. Although not a definite conclusion, 


this interpretation does nonetheless, present an interesting 


speculation. 
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Recommendations 


We should like to make the following recommendations, 


Since hindsight is better than foresight. 


1. More data are required along the Suffield east 
profile. This line should te extendéd past Swift Current to 
obtain first arrival information. Furthermore, the shot 
holes at Suffield should be drilled into Cretaceous rocks since 
we have found that shooting in the glacial till results in 


Very poor records. 


2. The filter setting should never be lower than 16 cps. 
sharp filter bandpass often results.in ringing in,the.filter 
circuit. Only for very distant shot (over 500 kms) should 


bower filter settings .be ,used. 


3. If at all possible, some sort of communications 
should be available between the shot point and the recording 


trucks. Such communication would eliminate much wasted time. 


4, A reverse profile should be shot in the mountains 
to make a definitive interpretation of the crustal structure 
under the mountains possible. Further work should also 
be done along the extended line, even if it means using 


multiple shot holes along this line. 
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Appendix 1 


Summary of Mathematical Principles Used 


Tiere ore Liairee Seclhians CoO This appendix. The rirst 
deals with the\theory of refraction seismology, for both 
the case,of horizontal layering and that’ of dipping layers: 
Refléctlon seismologyeis omitted since there were no studies 
carried out which involved reflected data. Secondly, a 
statistical section has been developed, showing those 
techniques used in calculating most of the information 
presenved in this thesis. And finally, “since some spectral 
analyses were carried out, it was felt that there should be 
ac least & cursory glancevsat the theory of autocorrelation 


and power spectra derived therefrom. 


Where the method of presentation is not my own, I 
have given credit where credit is due. This is particularly 
true of the sections dealing with statistics (taken in 
greater part from Steinhart and Meyer (1961)) and the section 
on power spectral analysis (taken in greater part from 


Lee (1960) and from Ellis (1964) ). 
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Mathematical Theory 


section A. 


Refraction Along N Horizontal Layers. 


Case 1. Single Refractor. 


Fig. (A-1) 


Consider a ray path as shown in the diagram above. 


total time along the path ABCD is given by 


Now 
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the first layer. It is the intercept, of the time distance 


line, at zero distance. This term may be considered as 
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being composed of two parts, one associated with the shot 

point, and the other with the detector. It is convenient 
mn cos. 0 

to call — the delay time; hence the intercept 


time is the sum of the delay times under the shot and 


detector. 


Case 2. N Refracting Horizons. 
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It can be easily shown that the total time T is given by 
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If one has refraction data along only one direction, then 
equation (3) allows one to calculate the various thicknesses 
h. of each layer, knowing the various velocities V5 of 


each layer from the travel-time graph. Furthermore, if one 
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has data from a double-ended profile, where the intercept 
times and apparent velocities are such that little dip is 
indicated in any of the horizons, then equation (3) can 
again be used to calculate the layer thickness under each 
Shot point. For a double-ended profile where the intercept 
times and apparent velocities are Significantly different, 
then the equations for dipping layers must be used. A 


discussion of the dipping layer case follows. 


it. Refraction Along N Dipping Layers. 


Case 1. Single Refractor. 
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We can then make the following substitutions 
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apparent velocity. In order to measure the dip, we must 


nave control along both directions of the profile. If 
this is indeed the case, then the apparent velocities 


yield the following information: 
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The’ total “travel, time equation is often seen in 


the form 
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Case 2. Two Refracting Horizons. 
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It is perhaps best to approach this problem from a 


Perrerent polnt of views 


In this case let us think in terms 


of apparent velocities, Vag being the apparent velocity 


measured at A, and 58 being the apparent velocity at F, 


for the at refractor. 
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From Snell's Law we can obtain directly, 
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be relegated to the N layer case to be discussed next. 


Case 3. N Refracting Layers. 
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By analogy to the double refracting case we can write down 


the following relations: 


Ms 
—— = gin (9, + 6,) 
Va UP as: I 
V ‘4 
=—— = sin (0 6. ) 
ne 1,nA i | 
Q -1 | Ye ot | 
rude Se Vy e 8 L,nn| 
lv | 
ee ee 
2) MB sin iV, sin "1, 0B) 


vd 
Vues 
om m+ 1 ; . 
= Mc hee! 
QnA sin v_ Si m, nA 
* fc" 
V 
-1 m+ 1 ; : 
= = epee 
Oy NB or Mee 


(16) 


7 
Li 


is, 


: *) 


“wwob SJI%W (eo Sw sé@Bge soe sso a8 > bans -W 
rae: 


: ono tale wan * a 


lads cart ales ody ' 
bess, 


furthermore 
0 = 2 
m,n, A m-1,n, A a on es 
(20) 
] = n 
m,n,5 nev) B - vy Oe oa | 
Finall h j i 
yy, the criticadangle net dh is given by 
Q - HeolaB rT VaRot hn (21) 
n-l,n e 
Q - 2 
aie) Z 
Dra = Sanne Eeeamenee = geP of Gop of 
nth jayer (225) 
V 
be = oo ='true'velocity in neh 
n-1,n layer (23) 


Having obtained these data for each layer, it is then 
possible to calculate the vertical thickness of each 
layer under the two shot points. By analogy to the 


single layer case, we may write the travel time equation as 
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T) = travel time along path from A 
R = distance from A 


hn, = thickness of layer “2 under “A, 


Using the zero distance intercept time and reducing equation (24) 


to a more familiar form we obtain 
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The analogous equation at B is written 
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Dipping layer theory was used for all the calculations 


carried out on the data shown in this thesis. 
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section B, 


Statistical Analyses. 


The following discussion is taken from Steinhart 


and Meyer (1961). 


ft. Least Squares Line FiCcing. 
Consider a line whose equation is T = tot Ax 
calculated from N points (t, x, ). The least squares 


estimates of the parameters A and vs may be found from 


the following equations 
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This latter equation is somewhat different from that 
quoted by Steinhart and Meyer (1961) since it was found 


that this was a more useful form for computational purposes. 


Having computed these parameters, it was then 
possible to calculate the thickness of each of the layers, 
their dip, and their true velocity using the various 


equations described in section A. 


Ii. Estimate of thesUncertainty of the»Slope A. 


Since the slope is used in the calculations of 
the crustal model, it seems obvious that we must have some 
means of describing the uncertainty of the slope. We 
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The value of Sy may now be employed to estimate S 


as follows: 


a oa (31) 


vone Blenificance of S> is not particularly clear as it 
now stands. A probability statement as regards the limits 
of uncertainty is perhaps the best way of expressing errors, 


and this is the type of statement we shall use. 


It may be shown rigorously that if w is the 
wrue value or the slope then A is distributed normally 
with mean Ww and variance la Acro pincée we are dealing 
with a finite sample, rather than resorting to tables of 
normal distribution to obtain our probability estimate, we 
mast resort to a table of Student's a distribution. Thre 
is the case since if mars 4 is normally distributed, then 
— will be ene according to Student's yi 
te ren We can then write the limits of the true 


value of w as follows: 
Vie kD NE (32) 


The values used in our calculations were for a probability 
of 0.1. This implies that nine out of ten time yY will 


have a value in the range as calculated using equation (32). 
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Tit, Uncertainty ofthe Time Intercept. 


Since we have used the time intercept for cal- 
culation of crustal thickness, a statement must be made 
regarding the uncertainty in by » ovstnhart and Meyer (1961) 
use the orthogonal form for the, variotié lines. Thus 


= o + AX when written in the orthogonal form becomes 


Tat + A (xX =x) (33) 
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Equation (36) can be used directly to find the 
uncertainty in the zero distance intercept time by setting 


X =O. Thus the confidence limits for t can be written as 
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This*°completes the discussion of the statistical methods 
used in analyzing the travel-time graphs and also in 


postulating our crustal model. 


Section C. 


The Measurement of Power Spectra. 


I. Fourier Analysis. 


Let us consider f(t) as being a periodic, piece-wise 


continuous function, whose integral 
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it may be expanded in a Fourier Series.of the following form 
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We note that F(n) is simply the Fourier Transform of f(t). 


This is a complete function as information regarding 


amplitude and phase angles is preserved. 


F(n) is in general complex, and is called the 
complex spectrum of f(t). Because the harmonic order n 
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The expression 


is of interest, where rs and fy, have the same fundamental 
angular frequency and t is a lag time between (-~,~) and 
is independent of t. An important property of this function 
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1 (n) being the complex conjugate of Rin). 
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We call (54) the Correlation Theorem for periodic functions. 
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The power spectra illustrated and discussed in 
this thesis were computed usingithe autocorrelation theorem. 
muse Tor completeness we discuss the autocorrelation funerion 
for oerLodiG Tunetilonsaand then simply wrive the analogous 


equations for random and quasi-random functions. 
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The left hand side of (55) is called the auto-correlation 


function and is defined by the symbol 1, (7). Therefore 


T,/e2 
04, (T) = 7 fit) fe + vj) dt (56) 
-T,/2 
Let ®, , (n) be the Power Spectrum in (55), 1.e:, 


eS A 


ot besaswelb bas bedattarhtl inion ale me Trea 7 : 
irmtioons noftelertosdJus dee a 15 
mlvoout motvalemosotus og aavowdb soomrguonel: yi 
ayogoisns adi efiaw ylqnta nent bas Batson ~ 


.enoiisonyt mobest-teasp bow mobnet "tf 


9 | Ma, Seg ee th 
+ menoerts noktefseton ond mt 


vent . t= 
S pi 


7m wi ae ve 
(a2) a [May 218 2 86, Geet) a ae 
he SY aye ; 

) Man ai ip 4 ib i | 
coiiaierteg-otus eft belles at A@) te pur mee bse v1 
axgtayot? <(?). rr? foamy sit yd bentteb at bri 06 

alt i= ath 


(2) ae (Pr + 9) 2 C9) ae Pag 7 


e- Set o (ee) Ge eee we 8 i 4 


—s 


te) 


A 26 


The Autocorrelation Theorem for a periodic function is 
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Let us now extend these equations to cover the random 


and quasi-random time series. 


Weiner Theorem for Autocorrelation. 
Por thercase jor the periodic function f, (t) we 
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Since os ta) is a real and even function, 
consequently its transform must be a real, even TUnCTL OI: 


we can thus write the above equations (62) as 
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ris Lc appears tnhav ®,, (w) is the power density spectrum 
of fj(t). It may also be shown that it is real, even and 


non-negative. 


TV. Computational Formulae. 
Summarized below are the computational formulae 
used to calculate the autopower spectra. This’ section is 


taken from Ellis (1964). 
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The Nyquist frequency of our data is SU ecpes . 
This frequency is sufficiently high such that the problem of 
aliasing is minimal in our Spectra, for the high-cut filter 
setting of the band pass was never higher than 16 cps and line 
balancing was always used to minimize any 60 eps components. 
see fig.( 2 ) for curve of the frequency response of the 


ViIF system for various filter SEvrings. 


Let L be a correlation tag index, Lb = G,- 41, -22, 
= (n-1) end tet J “be a frequency inde» forthe 
(m+ 1) equally spaced spectral estimates in the range 


OK Cams Wa ly - the following ouantities are computed. 


de AuGocorrelation 


n-L 
1 — 
ACL) er _ Xe ae ’ A(-L) = A(L) (63) 
ii 
[ue 20y Sle ee) » (1) 
ec. Autopower density spectra J Oe aa oe ve Segue ite 
(n-1) 
LJ 
Cae a D,(L) A(L) cos (==") (64) 
L=-(n-1) 


where D,(L) is a weighting function, called the Hanning 
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(Blackman and Tukey, 1958). This has the effect of 
convolving the Fourier transform of D, (t) LQ, (f) ~ the 
spectral window] with the raw Spectra. A simple and 
widely used lag window (whose use is called "hanning") 


is g.ven tn continuous form by 
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where T= maximum lag calculated. D,(t) and the 

corresponding Qo(f-f,) are shown in Fig. (A-6),- This 
window concentrates the main lobe of Qo(f-f,) near f. 
end tas rapidly falling sides. in digital calculation 


eee = MAG. 


Of A h 
a Ys 


7 7 

to goslts odd asd efAT he a mn oe LE 
oft - (7) 0) (7) @ %o wroledetd ners eit att i ovata 
bas ofgmte A .sitoesge nant ‘atts ddtw (Wosnt w “ietionae ef 
"“oitinsdt” befiso ef sae aod) wor tw gal been <a | 


vd sie? ‘@ound $rido nt ore 


( 


-< 


[ e>iet. Geary ge mg oo” 
(aay) Ay 

| T <4? O= mee . 

" i 


ed? bas (7), of heath a asl nwmtxan = 
atrit .(d-A) .aif at wore ous ( car ae 
1 aeon (3-1) @ to sdof them ent eedarda: mt 09 wobntw 

o ate oe cane 


notvsivolss ietiathb al .eebke gutted ir bias 


vee al 


Hanning lag window 


a: . OQ - a, 
te: 
on O45 
a 
O 


Hanning spectral window 


Change in 


4 
ae Ee da scale 
os, — Ove 
G 
0 O 
-0.02 
O a} 2 3 4 - 
ou, Lie se) 
Figure A-6. The lag and spectral window for the hanning 


estimate (after Blackman and Tukey, 1958). 
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